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*S . Abstract 

Lie group method provides an efficient tool to solve nonlinear partial differential equations. 
This paper suggests a fractional Lie group method for fractional partial differential equations. 
A time-fractional Burgers equation is used as an example to illustrate the effectiveness of the 
Lie group method and some classes of exact solutions are obtained. 
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1 Introduction 

Many methods of mathematical physics have been developed to solve differential equations, 
among which Lie group method is an efficient approach to derive the exact solution of nonlinear 
partial differential equations. 

Since Sophus Lie's group analysis work more than 100 years ago, Lie group theory has 
become more and more pervasive in its influence on other mathematical disciplines [1, 2]. There 
are, however, there few applications of Lie method in fractional calculus. Then a question may 
naturally arise: is there a fractional Lie group method for fractional differential equations? 

Some researchers investigated Lie group method for fractional differential equations in sense 
of the Caputo derivative and derived scaling transformation and similarity solutions [3-5] . Con- 
sidering the classical Lie group method, method of characteristic is used to solve symmetry 
equations. Recently, with the modified Riemann-Liouville derivative [6-8], we first propose a 
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more generalized fractional characteristic method [9] than Jumarie's Lagrange method [10]. Us- 
ing our fractional characteristic method, the generalized symmetry equations generating by the 
prolongation technique can be solved, and a fractional Lie group method was presented for an 
anomalous diffusion equation [9]. 

In this study, we investigate a simplified version of the fractional Burgers equations [5] 

u[ a) = Uxx + ul, x € (0, oo), 0<t, 0<a < 1, (1) 

and derive its group classifications. In order to investigate the local behaviors of the above 
equation, the fractional derivative is in the sense of the modified Riemann-Liouville [6-8] . 



2 Fractional Calculus and Some Properties 

From the viewpoint of Brown motion, Jumarie proposed the modified Riemann-Liouville deriva- 
tive [6-81, 



1 rl n l' x 



(2) 



where the derivative on the right-hand side is the Riemann-Liouville fractional derivative and 
n £ Z + . 

(a) Fractional Taylor series 

Recently, Jumarie- Taylor series [11] was proposed 

Here f{x) is a /ca-differentiable function and k is an arbitrary positive integer. 
Taking k = 1, f(x) is a a-differentiable function. We can derive that 

« = (4) 
r(i + a) 

(b) Fractional Leibniz product law 

If we set D^u(x) and D°v(x) exist, we can readily find that 

D«{uv) = u { ^v + uv {a) . (5) 

The properties of Jumarie's derivative were summarized in [11]. The extension of Jumarie's 
fractional derivative and integral to variations approach by Almeida et al. [12, 13]. Fractional 
variational interactional method and Adomian decomposition method are proposed for fractional 
differential equations [14, 15]. 

(c) Integration with respect to (dx) a 
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1 f'X ~i f'X 

oizm = — - / ( X -o Q - 1 /m = p7— r n / /(o«r,o<a<i. (6) 

r(aj j r(a + 1) Jo 

(d) Generalized Newton-Leibniz Law 

Assume D"/(x) is an integrable function in the interval [0, a]. Obviously, 



and 

D a 



r(i + a) 7o 

(e) Some other useful properties 



— ^— [ a D«f(x)(dxr = f(a)-f(0),0<a<l, 
1 (1 + a) Jo 

^jT^/(0(^=/(x) -/(«), 

f mmr = m, o<«<i. (9) 

JO 



£>a /3 = r (! + ft) /3-a < /3 < 1 



(7) 
(8) 



/ (a) ([z(t)]) = J^ (a) (*)> 0<a<l, (10) 



(11) 



J {dxf =x P . (12) 

The above properties (a)-(d) can be found in Ref. [11]. We must point out that f(x) should 
be differentiable w.r.t x in Eq. (10), and x 13 is an a order function in Eq. (11). 



3 A Characteristic Method for Fractional Differential Equations 

It is well known that the method of characteristics has played a very important role in math- 
ematical physics. Preciously, the method of characteristics is used to solve the initial value 
problem for general first order. With the modified Riemann-Liouville derivative, Jumaire ever 
gave a Lagrange characteristic method [10], in which the time- fractional order equals to the 
space- fractional order. We present a more generalized fractional method of characteristics and 
use it to solve linear fractional partial equations. 
Consider the following first order equation, 

o( x, f )^ + 6(x, t )^|a = C (x, f ). (13) 

The goal of the method of characteristics is to change coordinates from (x, t) to a new coor- 
dinate system (xq, s) in which the partial differential equation becomes an ordinary differential 
equation along certain curves in the x — t plane. The curves are called the characteristic curves. 
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More generally, we consider to extend this method to linear space-time fractional differential 
equations 

a(x,t)^^+b(x,t)^^=c(x,t), <a,fi<l. (14) 
With the fractional Taylor's series in two variables [11] 

similarly, we derive the generalized characteristic curves 

^ = c(x,t), (16) 

-a(x,t), (17) 



T(l + p)ds 
(dt) a 



. =b(x,t). (18) 
Eqs. (16)-(18) can be reduced as Jumaire's Lagrange method of characteristic if a = ft in [10]. 

4 A Fractional Lie Group Method 

In the classical Lie method for partial differential equations, the one-parameter Lie group of 
transformations in (x, t, u) is given by 

x = x + e£(x, t, u) + 0(e 2 ), 
i = t + er(x,t,u) + 0(e 2 ), 
u = u + e<j){x, t, u) + 0(e 2 ), 

where e is the group parameter. 

Use the set of fractional vector fields instead of the one of integer order 

V = £(x,t,u)D% + T{x,t,u)D? + (j)(x,t,u)D u , 0<a<l, < /3 < 1. (19) 

For the fractional second order prolongation PA 2 ^V of the infinitesimal generators, we 
proposed [9] 

— ~ v + ^ + ^ + ^ + ^ u + ^^Bf u ™ 

As a result, we can have 

Pr^V(A[u]) = 0, (21) 

onA[u] = 0. 
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In the time-fractional Burgers equation, Eq. (1), we only need to consider the case of the 
fractional order of space (3 = 1. Thus, the corresponding Lie algebra of infinitesimal symmetries 
is the set of fractional vector fields in the form 



V = £(x, t, u)D x + t(x, t, u)D? + <f>(x, t, u)D u . (22) 
We assume the one-parameter Lie group of transformations in (x, t, u) given by 

x = x + e£(x, t, u) + 0(e), 

u = u + ecj)(x, t, u) + 0(e), 

where e is the group parameter. 

The generalized second prolongation satisfies 



<D?u ^ v dD x u ^ V 

Using the following condition 



oDfu oD x u oDf a u ou xx oDfu x 



Pr (2) V(A[u]) = 0, A [it] = 0, (25) 

we can have 

(^-^ x -2«x^)| A[tt]=0 = 0. (26) 
The generalized prolongation vector fields are reduced as 

0* = £>? - (tffOArU " {D?T)Dfu, 

4> x = D x cf) - {D x £)D x u - (D x t )Dfu, (27) 
<T = Dl4> - 2(D x 0D 2 x u - [DlC)D x u - 2(D x t)D x D?u - (D 2 x t)D?u. 

Substituting Eq. (27) into Eq. (26) and setting the coefficients of u x u x a ^ , u x "\ u xx u x , u x 
and 1 to zero. Solve the equations with maple software, we can have 

£(x,t,u) = ci + c 4 x + 2c 5 r( *° a) + 4c 6 r ff a) , 

r(x, t, u) = c 2 + 2c 4 r(^y + 4c 6 r*(i+a) ' ( 2 8) 
4>(x, t, u) = c 3 - c 5 x + r 2 ( ^°) - c 6 x 2 + a(x, t)e u , 

where k[°^ = k xx . 

The Lie algebra of infinitesimal symmetries of Eq. (1) is spanned by the vector field 

\r d t/ _ d a T/ _ 9 t/" „, c? I 2t a d a 

^5 = -^ + r^5y&. ( 29 ) 

VQ ~ r(l+a) 3i r^(l+o) V x r(l+a)/3u' 
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and the infinite-dimensional subalgebra 



V k = k(x,t)e~ 



du' 



It is easy to check the two vector fields {V±, V 2 , V 3 , V4, V5, Vq} are closed under the Lie bracket 
[a, b] = ab — ba. In fact, we have 

[V h Vi] = (i = 0, 6), [Vi, V 2 ] = [V u V 3 ] = 0, [Vi, V 4 ] = -Vi, [Vi, V 5 ] = v 3 , 
[Vi, V 6 ] = -2V5, [V 2 , F3] = 0, [V 2 , V 4 ] = -2V 2 , [V 2 , V 5 ] = -2V U [V 2 , V 6 ] = 2V 3 - 4V 4 , 
[V 3 , V 4 ] = [V 3 , V 5 ] = [V 3 , V 6 ] = 0, [V 4 , V 5 ] = -V 5 , [V 4 , V 6 ] = -2V 6 , [V 5 , V 6 ] = 0. 
[V,V k ] = -V kx , [V 2 ,V k ] = -V kt , [V 3 ,V k ] = -V k , [V 4 ,V k ] = -v k ,, [V 5 ,V k ] = -V k n, 



[v 6 ,v k 



-V 



fc'"' 

2t° 



w V,p, p t,' _ T U ,_Z£_ k (a) ." _ 2t a h , r tflnH k'" - 4arf " h I , 4 ^" h [a) I ( r 2 I 2t " )k 
wneie/t — Xh, x -\- p( 1+a ) H > K — -p(l+a) x ^ ~ r(i+a) x ^W(T+a) t ^r(i+a)' 

Take the characteristic equation V5 as an example. The characteristic curve of V5 can be 



» 



given 



du 
de 

dx 2t a 



x, 



(30) 
(31) 

r(l + a)de ~ (32) 
Solve the above ordinary equations with the initial value u = u(x, t, e)| £ =o, x = x(e)\ £= o and 
t = t(e)\ £= Q. The one-parameter group Gi generated by the Vi(i = 1, 6, a) are given as 



de r(l + a)' 
{dtr 



0. 



5i : (z, 



r. u 



r(i+a) 

54: (£, r( - 1+Q ) 
55: (a; 



r(i+a) • 



56: (x, r(1+Q ), 

ffa: l x > r(l+a) ' 



( t a \ 

— > (X + £, p( 1+a ) j n J) 

~~ ^ V 2 '' rTl+a) £ ' U /' 

~~ ^ r(l+a) ' n + e J' 

/ £ t a e 2e \ 

v xe ' rTT+a) ' u /' 
— ?■ ^x -r r(i+ a ) ' r(i+a) > 

x t a 



24-a 



(33) 



l-4e. 



( x > r(TTa) ' log(e" + e/c)). 



r(l+a) 

w — 



xe 



r(i+a) " "°r(i+c«) 



l-4e 



r(i+a) 



Take a = 1 in the above classifications. We can derive the results for the case of integer 
order. Since gi is a symmetry, if u = f(x, r (i+ a ) ) ^ s a s °l u ti° n °f Eq. (1), then the following 
are also the solutions of Eq. (1) 
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u i — f{x e, r ( 1+Q ,)), 

u ? = f(x, nfr^y 

U3 = f(x, r(TT^y) + e ' 

u 4 = f(xe' £ , i^gy), (34) 

t( 2et OL t a \ i e 2 t Q1 

«5 — / - r(l+a) ' r(l+o) ^ + ~~ X£ ' 

Now we consider the applications of the above transformations. From u\ to U4, we can only 
obtain trivial solutions. Therefore, we start from the use of U5 

t a . . 2et a t a e 2 t a 

95- {x, — — — -, u)->(x + — — — t,^— — 7,^-777— — T-xe), (35) 
r(l + a) 1(1 + a) r(l + a) T(l + a) 

Assume 115^ = U5$(x, ^{l+a) ) = f(x, ^{l+a) ) is one exact solution of Eq. (1). Take 115^0 = c, 
where c is a arbitrary constant and also a trivial solution. We can get a new nontrivial exact 
solution as 

e 2 t a 

u 5,i = c + =— — r - xe. (36) 
T(l + a) 

Further more, continue this iteration process, we can derive a new exact solution of Eq. (1) 

2et a t a s 2 t a Ae 2 t a 

^ = U5 ' l(X " WTa) ' WTaj ] + WTa) - X£ = c ~ 2x£+ TOT*)' (37) 
Similarly, take «6,o = c, then we can have 



x 2 e , r t a 



1*6,1 = c-— — p l°g W 1 + 4e =^r— ■ — r, (38) 



and 



X 2 £ , I t a 



^6,2 = ^ 6 ,i( ta , — — ta )- t— P ~ log \/ 1 + 4£ rn wi (39) 

l + 4e r(1+Q) l + 4£p^y l+4£p^y y Hl + aj 



2ex 2 , I 8et c 



We can readily verify that 1*5^, 2x5,2, it6,i arid ^6,2 are four exact solutions of Eq. (1). 
Assuming the fractional order a = 1, the exact solutions we give here can be reduced as the 
exact iteration solution in Ref. [16] if we set the coefficients a = b = 1. 
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5 Conclusions 



Fractional differential equations have caught considerable attention due to their various appli- 
cations in real physical problems. However, there is no a systematic method to derive the exact 
solutions. Now, the problem is partly solved in this paper. The presented paper can be applied 
to other fractional partial differential equations and investigate their non-smooth properties. 
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